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EFFECT OF SUBCOOLING AND RADIATION ON FILM-BOILING 
HEAT TRANSFER FROM A FLAT PLATE 
by Thomas D. Hamill and Kenneth J. Baumeister 
Lewis Research Center 


SUMMARY 


A theoretical analysis of film boiling from a horizontal plate with subcooling and radi- 
ation is presented. The analysis is based on the postulate that the rate of entropy produc- 
tion is maximized. The general solution enables the overall heat-transfer coefficient 
h^t to be calculated when the following heat-transfer coefficients and subcooling group 
are specified: (1) saturated film-boiling coefficient, h^; (2) radiation coefficient, h ra( j; 
(3) turbulent free-convection coefficient for subcooled liquid, h tcl ; (4) subcooling param- 
eter, 9 = ( T - T, )/(T - T ), where T , T, , and T are the saturation temperature 

’ ' s o' ' w s” s’ b’ w 

of the liquid, the temperature of the bulk liquid, and the wall temperature, respectively. 
The general solution consists of a unique relation between (h^ - h rad>^Ib and 

(, 'racl " h tcl^/ h fb‘ 

A simple formula, which is accurate when | (h rac j - ^tcl^^fb' ^ was derived 
and is given by 


h tot " h fb + °* 88h rad + °* 12h tcl 6 


where the total heat flux q, . into the system is given by 

q tot = h tot( T w " T s) 


Theory predicts that film boiling is impossible for values of the group 
(h ra d " ht c i#)/hfb < -1. 27, that is, for strong subcooling. 


INTRODUCTION 


Radiation and subcooling are fundamental variables in film-boiling heat transfer that 



are particularly important in the quenching of metals. In quenching, radiation is a sig- 
nificant heat- transfer mechanism because the metal is initially at a high temperature. 

The temperature of the quenching bath is also important, because the properties of a 
metal resulting from an ice water quench can be markedly different from those obtained 
with a saturated quench. Subcooling, therefore, is of practical significance in quenching 
techniques. The purpose of this report is to extend the analysis of reference 1, in which 
film boiling of saturated liquids was considered with radiation neglected. The present 
analysis, which employs some of the concepts and results of reference 1, takes into 
account both radiation and bulk subcooling. It is hoped that the results will be of interest 
from a fundamental and a practical viewpoint. 

A concrete example is useful to clarify the general problem. A box with insulated 
sides, as shown in figure 1, is considered. The bottom of the box is made of steel plate 
which has been heated to a red glow. Cold water is then poured into the box. The plate 
is continually heated, keeping the temperature very high. The liquid that hits the plate 
is vaporized quickly, thus establishing an insulating layer of steam between the plate 
and the liquid. This condition is referred to as film boiling. The liquid that is evaporated 

is assumed to be simultaneously con- 
densed by a reflux condenser inserted 
into the vapor space above the liquid. 

This condenser cools the condensate 
below the saturation temperature before 
the condensate falls back into the liquid 
pool. 

The coolant temperature and flow 
rate in the condenser determine the 
various degrees of subcooling with re- 
spect to the saturation temperature that 
can be achieved. As a particular ex- 
ample, the wall may be at 1000° F, the 
interface at 212° F, and the bulk of the 
liquid at 50° F. The important point is 
that there are three temperatures that 
characterize the system: 

(1) Wall temperature, T w 

(2) Saturation temperature, T 

s 

(3) Bulk temperature of the liquid, 


Figure i. - Film boiling of water with bulk subcooling. Heat is transferred across the vapor 
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layer to the interface by conduction, convection, and radiation. Some of the heat that 
reaches the interface evaporates the liquid. The rest escapes into the bulk liquid by 
virtue of the temperature difference between the interface and the bulk liquid. Because 
of the agitation induced by the steam bubbles, heat transfer from the liquid- vapor inter- 
face is by turbulent free convection. 

In this report are considered two common problems that arise in discussions of film- 
boiling heat transfer: 

(1) The effects of subcooling on the results 

(2) When radiation is important 

From a design viewpoint, a method or formula is needed for computing the overall heat- 
transfer coefficient in those cases where both subcooling and radiation are present. This 
analysis indicates how the three heat-transfer coefficients (the saturated film-boiling 
coefficient, the radiation coefficient, and the liquid free-convection coefficient) are to be 
combined to yield the overall heat-transfer coefficient. The analysis is based on the 
postulate that the liquid-vapor interface attains an average configuration that maximizes 

the rate of entropy production of the system 
and surroundings. 


BASIC MODEL ASSUMPTIONS AND 
GOVERNING EQUATIONS 

An idealized model of film boiling is 
depicted in figure 2. A vapor film covers 
the entire plate with vapor domes spaced 
symmetrically. When a liquid is supported 
by a layer of vapor, as in film boiling, the 
liquid-vapor interface is inherently unstable 
in a gravitational field. At certain loca- 
tions on the liquid-vapor interface, vapor 
will break through and escape, under the 
influence of gravity, into the bulk liquid. 
Photographs of film boiling show that these 
escape points are dome-shaped cavities 
arranged in a cell-type pattern (see fig. 2). 

The major portion of the heat transfer 
to the liquid (see fig. 2) occurs across the 
thin portions of the vapor film. The vapor 



(a) Top view. Symmetrical distribution of cells on surface. 



(b) Side view. Direction of vapor flow into domes; axial 
coordinate, z; radial coordinate, r ; cell radius, Rp dome 
radius, R 1? gap thickness, 6. 

Figure 2. - Idealized model of film boiling on horizontal surface. 
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domes are so thick that essentially no heat is conducted into them. Physically, their 
function is to act as hydrodynamic sinks in which the generated vapor is collected. The 
assumed symmetrical distribution of these sinks implies that the velocity field in the 
thin film is radially symmetric. The proposed model is based on the assumption that 
there is some time-averaged or ensemble-averaged configuration of the system where all 
velocity, pressure, and temperature fields are at steady state, and that this statistically 
idealized configuration represents the average behavior of the actual system. 

The outer boundary of each unit cell in figure 2 is a circle of radius Rq; in reality, 
however, no plane surface can be mapped completely by nonoverlapping circles. A 
hexagonal boundary would be more accurate but this pattern leads to intractable mathe- 
matics. Mapping by circles does not account for a large number of semitriangular 
patches. The heat transfer to these curved triangular patches is assessed in an 
a posteriori manner as used in reference 1. The correction to the heat-transfer coeffi- 
cient which accounts for these patches is of the order of 10 percent. 

The heat-transfer coefficient can be obtained by solving the momentum and thermal 

energy equations for flow and heat transport in the thin annular vapor film of a single 

cell. The liquid- vapor interface is at the saturation temperature T , and the bulk liquid 

s 

is at a colder temperature T^. The bulk liquid is, therefore, subcooled with respect to 
the saturation temperature. The wall temperature T w is held constant. Radiation is 
assumed to occur across the vapor gap with negligible absorption within the vapor, and 
with no net radiative heat transfer within the liquid phase. All radiation is either absorbed 
or reflected at the liquid- vapor surface. The radiation rate depends on a "view factor" 
determined by the geometric shape of the liquid-vapor interface. An emitted photon can 
hit either the flat portion of the interface or the dome region. Reflection from these two 
surfaces would be different. A rigorous treatment of this problem would obviously be 
complex, and for simplicity, the view factor is assumed to be 1. This assumption com- 
prehends the first-order effects of radiation. 

Computation of the rate of heat transfer into the system for fixed wall, bulk, and 
saturation temperatures requires a solution of the momentum, continuity, and energy 
equations for all the boundary conditions. These equations are simplified by the following 
assumptions: 

(1) The inertia terms in the momentum equations are negligible, as justified in ref- 
erences 1 to 3. 

(2) The vapor is incompressible. 

(3) Convective effects in the vapor boundary layer are negligible. 

In addition, viscous dissipation is neglected, the flow field is assumed to be radially 
symmetric, and the physical properties are regarded as constant, even though they are 
evaluated at the film temperature T^ = (T w + T g )/2. With the previous assumptions, the 
governing differential equations in cylindrical coordinates are as follows: 
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Momentum: 


0 = 




0 = 




+ 1 + 9 w \ 

r3r 3z / 


Continuity: 


Energy: 


i-L (ru) + — = 0 
r 3r 9z 



Entropy production: 



= maximum 


( 1 ) 

( 2 ) 


(3) 


(4) 


(5) 


Equation (5) is the maximization postulate to be discussed more fully in the 
DISCUSSION OF MAXIMUM ENTROPY RATE POSTULATE section. The subscript 
’’universe" is used in a formal sense. In this report, the word universe means system 
and surroundings. The system is composed of those regions where temperature gradients 
exist, as well as any bubbles rising through the liquid. Thus, both the vapor layer and 
the liquid thermal boundary layer are part of the system. The surroundings consist of 
a heat source maintained at T w and a heat sink (the bulk of the liquid as well as a 
condenser) maintained at T^. The energy of the universe (system and surroundings) is 
constant, and thus, it can be considered an isolated system. In classical thermodynamics, 
the entropy principle applies only to isolated systems; that is, the entropy of an isolated 
system attains a maximum value in the final equilibrium state. 

In the mean statistical steady state, the total entropy of the vapor layer and the liquid 


thermal boundary layer is constant. That is, all the entropy generated in the vapor and 
liquid boundary layers per unit time is convected into the heat sink. Thus, in the steady 
state, only the reservoirs (source and sink) experience a net change in entropy with time. 
The entropy production of the universe is then given by 

M = o A 

' 'universe \ b 

and since T fe and T w are fixed temperatures, the entropy production is maximized if 
the rate of heat transfer is maximized. Therefore, the entropy production postulate 

implies that the rate of heat transfer will be a maximum subject to the following boundary 
conditions (see fig. 2, p. 3): 



z = 0 

u = 0 

w = 0 T = T w 

( 7 ) 

z = 6 

u = 0 

w = w. T = T_ 
0 s 

(8) 

Z = °° 


T = T b 

( 9 ) 

r = R 0 

u = 0 


(10) 


where Wg is the evaporation velocity at the interface. 

The momentum and energy equations are coupled at the liquid interface by a heat 
balance over the annular area: 

' PXW 6 = - k (f ) ■ h td (T s ' T b> + h rad< T w ‘ T s> <“> 

The term on the left side is the rate of release of latent heat (Wg is a negative quantity). 
The first term on the right is the rate of heat transport by conduction through the vapor 
film to the interface. The second term on the right is the heat escaping from the inter- 
face into the bulk liquid because of subcooling. The symbol h^ is the turbulent free- 
convection coefficient for heat transport within the liquid. The third term on the right 
represents the heat reaching the interface by radiation. The symbol h ra£ j is the radia- 
tion coefficient, a function of the plate and liquid emissivities and the wall and saturation 
temperatures. As stated previously, h ra£ j is assumed to be independent of the shape of 
the liquid -vapor interface. 
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Additional boundary conditions apply to the pressure field. A static-force balance on 
the annular liquid region requires that the average pressure on the flat portion of the inter- 
face be equal to the saturation pressure P : 

s 


2ir 



R r 


P(r, 5)r dr = tt(r 2 - R^P e 


( 12 ) 


This requirement must be met in order for the vapor to support the liquid. One further 
condition on the pressure field is that 


P(R r 6) = P g 


R 1 (p i - P)g 
Sc 



(13) 


This equation, which is derived in reference 4, relates the pressure at the entrance to 
the dome P(Rp 6) to the saturation pressure and to the difference in head within the dome 
corrected for surface tension effects caused by curvature. The pressure at the radial 
entrance must always be less than the system pressure in order for the dome to function 
as a hydrodynamic sink. A minimum dome size can be computed from equation (13) by 
letting P(R.,6) = P . This requires that 

-L b 


R i> 



2g g c 

i ~ p)s 


(14) 


If Rj was less than this value, the pressure inside the dome would be greater than the 
pressure within the annular film and flow into the dome could not occur. 


GENERAL METHOD OF SOLUTION 
Momentum Equations 

The momentum equations that apply in this analysis are solved in reference 1. The 
essential result is that in order to support the liquid phase at a height 6 above the plate, 
the evaporation velocity Wg required is given by 
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( 15 ) 


w g = -l/3 2 5 3 
0 2 


where P , a function of geometry, gravity, and physical properties, is 


p = 


g c 

S^Rq 


° 4 gp 2 f 


1 - f 


+ (^ + 


In f 


(16) 


and f, the ratio of the dome radius to the cell radius, is 



(17) 


Equation (15) is a functional relation between the evaporation velocity and the thick- 
ness of the gap that follows solely from momentum considerations. No matter how the 
vapor is generated, equation (15) must always apply. For example, a relation identical 
in form to equation (15) is derived by Whitney (ref. 5) and verified experimentally by 
Pearson and Bradfield (ref. 6) in their study of the gas-levitated disk. In the gas- 
levitated disk experiment, a stream of air is blown through a hole in a small solid disk 
at a rate high enough so that the disk floats on a cushion of air. The relation between the 
velocity of air flow and the distance the disk equilibrates above the plate is given by an 
equation identical in form to equation (15). While in the levitated disk problem the vapor 
flow is generated mechanically from a source of compressed air, in film boiling the 
source of vapor is evaporation at the interface due to the heat transfer across the gap. 


Energy Equations 

To keep the mathematics simple, the effect of heat capacity Cp (or equivalently of 
convective effects in the vapor layer) has been ignored for the present. This simplifies 
the solution of the energy equation to a linear drop of temperature across the gap. This 
defect is ultimately corrected by attributing the effect of Cp to the saturated film- 
boiling coefficient derived in reference 1. Thus, the temperature profile in the annular 
vapor layer is given by 
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(18) 


T <*> = T w - < T w - T s> j 


Substituting equations (18) and (15) into the interface energy balance (eq. (11)) yields 


i pwV = 


T s> 


h tcl< T s - T b> + h 


(T - T 1 
rad' w s' 


(19) 


Equation (19) is a quartic equation for the gap thickness in terms of the system param- 
eters. 

It is emphasized that equation (19) is a heat balance over the annular area only. The 
total rate of heat transfer from the plate per unit cell is given by 

Qtot = * R o| < T w - V* 1 - f2 > + h rad < T w ' T s>] + ^dome < 20 > 

In equation (20), the heat flux at the wall is being considered. In appendix B, the heat 
flux by conduction under the dome m is shown to be negligible in comparison to con- 
duction under the annulus if the gap thickness is much smaller than the dome radius. 

Thus, if conduction under the dome is negligible, the total heat-transfer coefficient per 
unit cell is 

h tot = = | (1 - f 2 ) + h rad (21) 

’ R 0< T w- T s> 


Optimization Procedure 


The maximization postulate requires that the amount of heat transferred compatible 
with the boundary conditions be a maximum; that is, that and Rq adjust themselves 
in such a manner that h tQ ^. is maximized: 



( 22 ) 
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(23) 



= 0 


From a mathematical viewpoint, the problem is now completely formulated. The depend- 

o 

ent variables (unknowns) at this point are /3 , 6, h tot , R Q , and Rp and there are five 
independent equations (eqs. (16), (19), and (21) to (23)). The maximization postulate, in 
essence, specifies what values of Rq and R^ will exist physically. Without the postu- 
late, empirical equations (or some other theory, such as hydrodynamic stability consid- 
erations) would have to be used for Rq and Rp From this point on, the solution of the 
equations is merely a matter of algebraic manipulation and numerical techniques. Details 
of the algebraic manipulations and rationale of the numerical procedure are shown in 
appendix C. 


Dimensionless Groups and Parametric Solution 

The parameters of the problem are combined into dimensionless groups. A dimen- 
sionless gap thickness r] is defined as 


(24) 


where 


/ cr Sr 

l =+ — V ( 25 ) 

Y (p z - p)s 

The characteristic length l emerges naturally from the analysis presented in appendix C. 
The optimal dome and cell radii are multiples of l. 

Radiation and turbulent free-convection Nusselt numbers are defined, respectively, 

by 


N 


W 


rad 


N. 




tel 


(26) 

(27) 
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where k is the thermal conductivity of the vapor in all Nusselt numbers in this report. 
The film-boiling Nusselt number in the absence of subcooling and radiation was 

derived from the application of the maximization postulate in reference 1 and is given by 


N fb = 


h fb* 


0.41 Ra 


a/4 


C P (T W - T s> 


(28) 


where Ra is the Rayleigh number and X* a modified latent heat; both are defined 
further in appendix A. In addition, it is convenient to define the subcooling parameter as 


e = 



(29) 


The symbol 6 is the ratio of liquid subcooling to vapor superheat. For a particular prob- 
lem involving subcooling and radiation, these three Nusselt numbers and the subcooling 
group can be computed a priori. From a design viewpoint, a method for relating these 
four groups to an overall Nusselt number for the system is required. The overall Nusselt 
number N^ is defined by 


N. 


tot 




(30) 


where 


q tot “ h tot< T w " T s> 


(31) 


As shown in appendix C, the general solution was reduced to two parametric equations 
(see eqs. (C35) and (C38)) of the form 


N tot ~ N rad 


N. 


= <P( f) 


fb 


N rad - N tcl 0 


N. 


= m 


fb 


(32) 


(33) 
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where cp and \f/ are complicated functions of f. If values of f are substituted into the 
right sides of equations (32) and (33), a numerical correspondence between the two left 
sides of the equations is generated. In other words, the group (N^ - N rad )/Nj b is a 
unique function of the group (N rad - which can be calculated from equations 

existing in the literature. Thus an overall Nusselt number can be obtained. 


RESULTS 


As stated in the preceding section, the general solution consists of a unique relation 
between the group (N tot - N rad )/ N fb and th « group (N rad - N tcl 0 )/N fb - The general 
trends obtained after numerically evaluating equations (32) and (33) by substituting values 
of f are given in table I. 

A graph of (N tot - N rad )/N fb 

plotted against (N rad - N tcl 0 )/ N fb 
is given in figure 3, and this graph- 
ical function is the general solution 
to the film-boiling problem with radi- 
ation and subcooling. Table III (in 
appendix C, p. 31) contains selected 
tabular values, including values of 
f, corresponding to the two groups. 

A peculiar anomaly, inherent in 
the parametric equations, arose in 
computing the general solution. When values of f less than about 0. 4 were substituted 
into the parametric equations, a new branch solution was generated (indicated by the 
dotted line in fig. 3). This branch solution obtained for values of f higher than 0. 4. 

Since this branch solution leads to a smaller heat transport, it was first thought to be 
physically meaningless; however, this branch turned out to be a pure free-convection 
line. That is, except for slight numerical errors introduced by the correction for the 
triangular paths, the equation for the dotted line is 


TABLE I. - HEAT-TRANSFER REGIME 


N rad N tcl* 
N fb 

Characteristic 

N tot - N rad 
N fb 

Positive 

Radiation 

<1 

Zero 

No subcooling or radiation 

1 


(exactly balanced) 


Negative 

Subcooling 

>1 

<-1.27 

No film boiling possible 

Undefined 


N tot " N rad 


N rad - N tcl 0 


N. 


fb 


N. 


fb 


(34) 


or 


N tot = N tcl 0 


(35) 
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Figure 3. - General graphical solution of effect of subcooling and radiation on 
film boiling. 


and the total heat flux is then 


or simply 


q tot h tot< T w " T s> h tcl( T w ” 


T s> 



q tot “ h tcl< T s ' T b^ 


(36) 


(37) 


What this means, then, is that the optimization procedure allows two solutions: one 
stable (film boiling), and one metastable (pure free convection). 

However, the analysis makes an even more remarkable prediction. The intersection 
of the film-boiling line with the free-convection line places a theoretical bound on the 
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possibility of film boiling in the presence of subcooling. That is, theory predicts that 
film boiling is impossible whenever the value of the group (N rad - ^tcl^^fb is ^ ess 
than - 1. 27 (i. e. , strong subcooling). For values of this group below - 1. 27, something 
must happen to film boiling; either it reverts to pure convection or to nucleate boiling. 

No physically meaningful value of f can be found for values of the group below -1. 27. 
However, it is apparent from examining figure 3 that the film-boiling line becomes tan- 
gent to the free-convection solution at - 1. 27. Therefore, for values of (h ra( j - h tcl 0 )/ h fb 
below this critical value, the free-convection line is the correct and only solution. 

At extremely high positive values of (h a( j - h.^J/h^ corresponding to an f value 
of 0. 793, the gap thickness becomes zero and the model no longer applies. As 
(\ad ” ^tcl^'Afb approaches infinity (in the vicinity of f = 0. 793), the group 
< h tot - h ra d)/hfb approaches zero. This implies that h^ = h ra( j in the limit of infinite 
radiation, which is physically correct. Why the limit of infinite radiation should corre- 
spond to an f value of 0. 793 is not known. 

The optimum dome size Rj, which is given by equation (C6) in appendix C, is inde- 
pendent of subcooling and radiation: 



6ag r 


] (pl ~ P)g 


(38) 


Theory predicts, therefore, an upper limit to the spacing of the vapor domes when 
f = 0. 4 which is given by 


D 


6ag, 


0, max o. 4 y (pj - p)g 


= 12. 24 1 


(39) 


For purposes of comparison, the most dangerous wavelength computed from hydrodynamic 
stability theory is given by 10. 85 1 (ref. 4). The two numbers are of the same order of 
magnitude. The vapor domes become more widely spaced with subcooling and more 
closely spaced with increased radiation. The average size of the vapor domes, however, 
should remain constant under all conditions. 

A simple formula for the overall heat-transfer coefficient was obtained by drawing a 
tangent to the graphical function at (N ra( j - N^^SJ/N^ equal to zero. The formula is 


h tot - h fb + °' 88h rad + °- 12h tcl 


( T s - T b> 
< T w - T s> 


(40) 
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This linear approximation ought to be valid for many practical conditions encountered. 

It is quite accurate in the range | (N ra( j - 0)/NfjJ < 0. 5; that is, if the film- boiling 
coefficient is twice the value of |h ra( j - h^j#], equation (40) is accurate. A numerical 
example of the film boiling of nitrogen is presented in appendix D to illustrate the theory. 

The simplified formula (eq. (40)) is similar to one derived by Bromley (ref. 7) for 
forced-convection film boiling from the outside of pipes. For fully developed turbulent 
(liquid) flow with no subcooling but with radiation, Bromley recommends the equation 

h tot = h fb + ~ h rad 

where h^ is the saturated film-boiling coefficient appropriate to the pipe geometry and 
flow conditions in the absence of radiation. 

The coefficients 0. 88 and 7/8 are so close as to lead to the speculation that equa- 
tion (40) may be a general relation for many geometries and flow conditions, if the cor- 
rect individual cpefficients are used. 


DISCUSSION OF MAXIMUM ENTROPY RATE POSTULATE 

The present analysis and its predecessor (ref. 1) rely heavily on the postulate that 
the system will attain a state in which the rate of entropy production of the universe is a 
maximum. The universe (system and surroundings) is always an isolated system, and 
it is to isolated systems that the postulate is limited in analogy with classical thermo- 
dynamics. Since the word "entropy" has been used, it is natural to look to equilibrium 
thermodynamics for clarification. 

The central problem of classical thermodynamics is to predict the equilibrium state 
of a system resulting from the removal of a barrier. For example, a rigid insulated 
cylinder separated into two parts by a barrier is considered. On one side is a gas at 
high pressure, on the other a gas at low pressure. The problem is to determine what 
final state will result for the overall system when the barrier is removed. 

Callen (ref. 8, p. 24), in discussing this problem, states: 

. . . the tentative postulation of the simplest formal solution of a problem 
is a conventional and frequently successful mode of procedure in theoretical 
physics. What then is the simplest criterion that can reasonably be imagined 
for the determination of the final equilibrium state? From our experience 
with many physical theories we might expect that the most economical form 
of the equilibrium criterion is in terms of an extremum principle. 


15 


I 


The extremum principle of thermodynamics is that the entropy of a fixed- energy system 
is a maximum in the final equilibrium state. 

In other words, the simplest formal solution to problems involving uncertainty or 
randomness is to maximize or minimize some quantity; the idea is that the extremum 
principle generates the required number of equations to make the problem determinate. 
Although there is no obvious connection between equilibrium and nonequilibrium processes, 
both would appear to be amenable to the same type of solution. Something is known about 
the problem; for example, in the first problem, perhaps the energy is constant, or in the 
film-boiling problem, perhaps the boundary temperatures are fixed. But other boundary 
conditions cannot be known in sufficient detail to make the problem determinate. It is 
logical to hope that a mathematical procedure that worked for the first problem will also 
work for the second. 

This kind of thinking led Malkus (refs. 9 and 10) to a possible breakthrough in the 
study of turbulence. First, Malkus considered turbulent natural convection from a hori- 
zontal surface (ref. 9). He postulated that the turbulent fluid motions would attain an 
extreme state in which the maximum amount of heat compatible with the boundary condi- 
tions would be transferred. By using this postulate, he deduced the relation of the Nusselt 
number against the Rayleigh number to the one-third power, an expression for the tem- 
perature profile, and an estimate of the mean- square velocity distribution - a significant 
achievement considering the a priori approach. Zuber (ref. 11) subsequently adapted the 
results of Malkus to the problem of nucleate boiling. The present analysis and reference 1 
apparently complete a series of papers on heat-transfer processes from horizontal sur- 
faces in which the heat-transfer rate is a maximum (see table II). 

Malkus (ref. 10) then used the same postulate (maximum rate of entropy production) 
to determine the turbulent-velocity profile in a flat channel. He deduced profiles similar 
to the defect laws of von Karman. Recently, Nihoul (ref. 12) applied the Malkus theory 
to magnetohydrodynamic turbulent channel flow. 

One of the differences between the film-boiling analyses and the Malkus papers is the 
degree of complicacy of the mathematical techniques. Malkus considers turbulent flow. 
Associated with this condition are stability arguments and complicated mathematics. In 
film boiling the vapor flow is laminar, so the mathematics are simple. In fact, the tur- 
bulence in film boiling exists only at the liquid-vapor interface and not in the vapor flow 
itself. Yet, the Malkus principle still seems to be correct: maximize the rate of entropy 
production to make the problem determinate. 

Just as Carnots’ reasoning on steam engines was the beginning of a much wider 
principle (Second Law), it may well be that the Malkus postulate has wide applicability in 
engineering problems characterized by macroscopic uncertainty. Turbulent flows are 
time dependent, yet average values emerge. Many other engineering problems, such as 
two-phase flows, are time dependent, and for these problems the maximization postulate 
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TABLE E. - PROBLEMS SOLVED BY MAXIMUM ENTROPY RATE POSTULATE 3 


Process 

Author 

Results 

Source 

Shear turbulence 

Malkus 

Universal velocity profile 

Ref. 10 

Magnetohydrodynami c 
(MDH) turbulence 

Nihoul 

Effect of magnetohydrodynamic variables 

Ref. 12 

Turbulent free 
convection 

Malkus 

Heat -transfer coefficient 
Temperature profile 
Eddy velocities 

Ref. 9 

Nucleate boiling 

Zuber 

Heat -transfer correlation 
using Malkus theory 

Ref. 11 

Film boiling: 
Flate plate 
Saturated 

Hamill and 
Baumeister 

Heat-transfer coefficient 
Optimum dome size 

Ref. 1 

Subcooling and 
radiation 


Optimum spacing of domes 
Effect of subcooling and radiation 

Present 

report 

Horizontal tube 


Effect of geometry 

Ref. 16 


21 

Physical problems involving rate processes that are governed by unknown or chaotic 
boundary conditions are resolved by maximizing the rate of entropy production of 
the universe. 


may work, whether or not the flow itself is turbulent. Time dependence or uncertainty 
seems to be the key condition. In any event, the maximum entropy rate postulate is a 
technique worth trying whenever an analysis is hindered by a lack of information. 

CONCLUDING REMARKS 

The problem of film boiling with subcooling and radiation has been treated by means 
of the maximum entropy rate postulate. The overall heat-transfer coefficient for the 
process can be computed by 


h tot " h rad + ^ h fb 

where cp is a unique function of the individual heat-transfer coefficients (radiation, h ^ 
turbulent liquid free convection, h. p saturated film boiling, h ffe ) and a subcooling pa- 
rameter 9 = (T - T,)/(T - T ). That is, cp is given functionally by 

b u w s 

- _ J h rad " h tcl 

V h fb 
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which is given graphically in figure 3 (p. 13). A simple formula that approximates the 
numerical results is given by 


h tot “ h fb + °* 88h rad + °* 12h tcl 


T - T u 
s b 

T - T 
w s 


which is valid for 


h rad ~ h tcl g 
h fb 


<0.5 


The analysis predicts that film boiling is physically impossible for strong subcooling 
in particular, when 


h rad - h tcl g < 2? 

^fb 

This prediction and the heat-transfer relations should be verified experimentally. 

The results obtained could possibly be valid for systems other than pool film boiling 
on a flat plate. For example, they may be valid for film boiling in other geometries if 
the appropriate individual heat-transfer coefficients are employed. The same thing may 
be true for forced convection with film boiling. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 15, 1966, 

129-01-09-04-22. 
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APPENDIX A 


SYMBOLS 


A 

B 


'P,l 


D 


0, max 


F(f) 


g 

Sc 

h fb 


‘rad 


‘tel 


tot 


1/RaA 

< N rad - N t cl # )/ RaA 

specific heat of vapor, 

Btu/lb mass 

specific heat of liquid, 

Btu/lb mass 

maximum predicted diameter 
of cell (eq. (39)), ft 

function of f (eq. (C8)) 

ratio of vapor-dome radius to 
cell radius, R^/Rq 

2 

local value of gravity, ft/sec 

Newton's law conversion factor, 
32. 174 (lb mass)(ft)/ 

(lb force) (sec 2 ) 

saturated film-boiling coeffi- 
cient, Btu/(sec)(ft 2 )(°R), 


0.41 


* 


k X pg(p t - p) 


1 1/4 


“< T W - T s>' J 


radiation coefficient, 

Btu/ (s ec) (ft 2 ) (°R) 

turbulent-free-convection 
coefficient for subcooled 
liquid, Btu/(sec)(ft 2 )(°R) 

overall heat-transfer coeffi- 
cient (eq. (31)), 

Btu/ (sec) (ft 2 ) (°R) 


N. 


fb 


N 


rad 


N. 


tel 


N. 


tot 


N 

P 

Pr, 


^dome 

^tot 


thermal conductivity of vapor, 
Btu/(sec)(ft)(°R) 

thermal conductivity of liquid, 
Btu/ (s ec) (ft) (°R) 


characteristic length, 
, ft 




IS (P t ~P) 

saturated film-boiling Nusselt 
number, h^Z /k 

radiation Nusselt number, 
h rad^ 

turbulent free- convection 

Nusselt number for subcooled 
liquid, h tcl ZA 

overall Nusselt number for 
system, h to) /A 

N tot “ N rad’ see < c29 ) 

2 

pressure, lb force/ft 

Prandtl number for liquid, 

C p, l M l 

saturation pressure, 
lb force/ft 2 

rate of heat transfer under 
dome, Btu/sec 

total rate of heat transfer into 
system, Btu/sec 
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^dome 

heat flux under dome, 
Btu/ (hr) (ft 2 ) (°R) 

z 

distance from plate to arbitrary 
point in vapor, ft 

^tot 

total heat-flux into system, 
Btu/ (hr) (ft 2 ) (°R) 


thermal expansion coefficient for 
liquid, (°R ) -1 

Ra 

Rayleigh number for film boiling, 
Z 3 gP C p(P* ~ P)/P k 

/3 2 

parameter given by eq. (C7), 
l/(ft 2 )(sec) 

R o 

R 1 

radius of cell, ft 
radius of vapor dome, ft 

6 

t(f) 

vapor-layer thickness, ft 
correction factor, eq. (C7), 
dimensionless gap thickness, 

R o 

optimal cell radius, Rj/f, ft 


6/1 


optimal dome radius, 
./ 6ag c 

e 

subcooling parameter, 
< T s - T b>/< T w - T s> 


yg(p* - p) ’ 

A 

enthalpy ratio, X*/C p (T w - T g ) 

r 

» 

radial coordinate (see fig. 2 ), 
ft 

X 

latent heat of vaporization, 
Btu/lb mass 

S 

entropy of universe, Btu/°R 

X* 

modified latent heat, 

T 

temperature, °R 

o 


xfi + M 0 p (T w- T s>, 
V 20 p X J 

T b 

temperature of bulk liquid, R 


Btu/lb mass 

T f 

film temperature, 

T f = (T w + T s )/2 > ° R 

M 

viscosity of vapor, 
lb mass/(ft)(sec) 

T s 

saturation temperature of liquid, 
°R 

H 

viscosity of liquid, 
lb mass/(ft)(sec) 

T w 

wall temperature, °R 

t 

perpendicular distance from wall 
to point on dome surface, ft 

t 

time, sec 

P 

3 

density of vapor, lb mass/ft 

u 

w 

radial velocity of vapor, ft/sec 
axial velocity of vapor, ft/sec 
evaporation velocity at the inter- 

Pi 

a 

O 

density of liquid, lb mass/ft 
surface tension, lb force/ft 

0 

face, ft/sec 

(P 

<?(« = (N tot - N rad )/N fb (see 

y 

dummy variable used in ap- 


also eq. (C38)) 


pendix B 


^ (f) = (N rad ‘ Ntcl^/Nfb (see 
also eq. (C35)) 
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APPENDIX B 


JUSTIFICATION FOR NEGLECTING CONDUCTION UNDER VAPOR DOME 

Conduction under the dome can be estimated by assuming a linear drop of temperature 
from the plate to any arbitrary position on the curved interface. Thus, the temperature 
profile at any radial position under the dome will be 


T(z) = T w - (T w 



(Bl) 


where £ is the perpendicular distance from the wall to any point on the surface of the 
dome (see fig. 4). The symbol £, which is a function of radial position, is given by 


€ = 6 



(B2) 


The total heat flow under the dome is obtained by integrating the heat flux over the dome 
area: 



This integral can be evaluated exactly by letting 


(B3) 



duction under dome. 



and noting that 
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or 


Therefore, 


dr 2 = 2r dr 




y" 2 dy 



= 27Tk(T w 

= 2tfk(T w 



- - 6 In y 

y 


1/6 

1/(5 + Rj) 


- T s> 


6 + R, - 6 - 6 In 


1 

6 > 

1 


6 + R^ 


= 27Tk(T w - 


T ) 
s 7 



6 In 


5 + V 


(B4) 


Substituting this value of Q dome for 7rR^q dome in equation (20) yields the total rate of 
heat transfer at the wall as 


Qtot - «0< T w - T Jj - ‘ 2 > + h rad] + 2l *< T w - T s> 


R 1 + 6 In 


6 


,6 + R, 


(B5) 


22 



I 


The overall heat- transfer coefficient h^ is then 


h tot = 


Q, 


2kR 1 


tot = k (1 _ f 2) + h -n + M m - 6 - 

ttr2(T w -T s ) 5 r5 Rn \ 5+R V 


(B6) 


v 0 0 


or 


U - k /! A , u , 2kf , 2kS / 5 

‘' t ot-< 1 -‘) + h rad + _ + _ln^_ 


(B7) 


l' 


It is seen that if Rq and R^ are much greater than 6, then 


h tot = 7 (1 - f2 > + h rad 


(B8) 


6 


which is equation (21) in the main text. A question might arise as to the order of mag- 
nitude of the third term in equation (B7) with respect to the first. Dividing one by the 
other yields 


2kf 

R 0 „ _ 6 _ 

^ (1 - f 2 ) R 0 

6 

-2 

Calculations indicate that 6 /R q is of the order of 10 ; therefore, the third term can 

be safely dropped. 
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APPENDIX C 


METHOD OF SOLUTION 
Algebraic Manipulations 

A simultaneous solution of equations (19), (21), (22), and (23) in conjunction with the 
o 

expression for /3 (eq. (16)) is required. The maximization with respect to Rq is con- 
sidered first. Equation (22) can be rewritten in terms of new variables as 



[*W R 1’ R 0>] 



(Cl) 


(see appendix E). From equation (21), the following equation can be deduced: 



= 0 


(C2) 


This condition requires that the gap thickness be minimized with respect to Rq in order 

to maximize the heat transfer. Equation (19) gives an implicit equation for 6 in terms 

2 

of Rq. In equation (19), |3 is a function of Rq as shown in equation (16). Taking the 
derivative on both sides of the interface energy balance (eq. (19)) with respect to Rq and 
solving for ^86 /3Rgj give 



V 2 3R o/ £ 

W T w- T s )- h tcl< T w- T s )- 2 P x ' 3263 


(C3) 


The only way for this equation to be satisfied for all arbitrary values of the parameters is 
for 



(C4) 
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Taking the derivative of /3 (eq. (16)) with respect to Rq, setting this equal to zero, and 
solving for R Q yield 


R* 1 / 

° f y (Pj ■ p)s 


(C5) 


or 


R 


* 

1 


6cg „ 


ip t - p)g 


(C6) 


The optimum vapor-dome radius Rj given by equation (C6) is identical to the result ob- 
tained in reference 1. The determination of the optimum f is not as straightforward - a 

O 

numerical approach is necessary. Substituting Rq into the expression for /3 (eq. (16)) 
yields 


p2 _ 2 ebi - p ) 

P* F(f) 


(C7) 


where Z is a characteristic length given by equation (25). The motivation for its intro- 
duction is based on Rj. A function only of f is given by 


F(f) = (-3)6 


3/2 


(i - <1 

4 


+ l^r + lnf 


f 2 (l - f 2 ) 


(C8) 


A characteristic length given by R ^/ ^6 emerged quite naturally from the first maxi- 
mization. This characteristic length is used to make the equations dimensionless. Mul- 
tiplying equation (21) by 1 / k gives the dimensionless equation 


N. 


1 - f‘ 


tot 


+ N 


rad 


(C 9) 


The interface energy balance (eq. (19)) becomes 
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4 


22 _ 

F(f) 


k(T - T )p 
' w s'^ 


pXg(p z - p)Z' 


+ V 


rad v w 


< T - - T *> - h tcl( T s - T h) 


pMp z ~ PW 

P 


(CIO) 


The characteristic Rayleigh number for film boiling is given by 


Ra = 


l gpC p (p t - p) 

pk 


Employing this dimensionless group yields 


t) 4 = F(f) 


> + i N -«- N 


RaA 


RaA 


where 


(Cll) 


(C12) 


6 = 


T - T . 
s b 

T - T 
w s 


is a subcooling parameter and replacing X by X gives 


A = 


C (T - T ) 
p' w s' 


(C13) 


The maximization with respect to f will now be undertaken. At this point, there are 
three unknowns (r\, f, and N^) and three equations (eqs. (C9) and (C12) and the maximi- 
zation with respect to f) to be considered. Thus, for maximum heat transport, 


dN tot = -(1 - f 2 ) dr? _ 2f _ 0 

df 2 df 7? 

V 

see equation (E3). Solving for dp/df yields 


(C14) 
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(C15) 


dr/ _ _ 
df 


2 f 


1 - f 


2 


V 


Before taking the derivative of equation (C12) with respect to f, it is convenient to let 

A = — — (C16) 

RaA 


and 


B = 


N rad - N tcl g 
RaA 


(C17) 


Thus, the energy interface balance becomes 

r, 4 = (A + B 7 ])F(f) (C18) 

From this point on, F = F(f). Taking the derivative of 77 with respect to f yields 


477 3 = (A + B 77 ) — + BF 

df df df 


Solving for drj/df gives 


d?7 _ 
df 


(A + Br,) ^ 
df 

477 3 - BF 


V 


4 dF 
df 


(4t7 3 - BF)F 


Equating (C15) to (C20) gives 


4 dF 

77 — 

2 f df 

77 = — 

1 - f 2 (4t7 3 - BF)F 


(C19) 


(C20) 


(C21) 


The goal of these manipulations is a set of parametric equations relating A and B to f. 
Solving for B from equation (C21) gives 
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(C22) 


B = 77 



1 - f 2 1 dF 
2f f 2 df 


but from equation (C18) 


A = 



Substituting equation (C22) into equation (C23) gives 


A = rj 


4 


3 _ 1 - f 2 - _1_ dF 

F 2f „2 df 
r t 


Thus, 



A 1 / 4 

1 - f 2 J_ dlA 1//4 

2f f 2 df J 


Cubing equation (C25) and substituting for 


3 




in equation (C22) give 


A 3 / 4 


_4_ [ 1 - f 2 1 dF 
F 2f p 2 df 



2 \ 3 / 4 
1 - i_ 1 dF | 

2f F 2df/ 


Substituting equation (C25) into the expression for the total Nusselt number (eq. 
gives 

/ 2 \ l/ 4 

- f 2 \ L 1 . 1 ~ f J_ dF \ 


(1 - n 


N. 


tot 


2f ^2 df / 

. 7 +N 


1/4 


rad 


(C23) 


(C24) 


(C25) 


(C26) 


(C9)) 


(C27) 
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Equations (C26) and (C27) constitute a pair of parametric equations for N^, the indi- 
vidual Nusselt numbers, and the value of f. That is, if an arbitrary value of f is sub- 
stituted into equations (C26) and (C27), unique values of B/A 2 ^ 4 and A^' 4 (N^ 0 £ - N ra( j) 
result. Numerically computing these values for the entire range of admissible f values 
results in a complete correspondence of to the individual Nusselt numbers. This 

constitutes the general solution. It is convenient to define 


N* = 1 " f 


1/4 


3 _ 

F 


1 - f 
2f 


\ 1/4 

J_ dF \ 

F 2 */ 


Therefore, 


N tot = N ’ 


+ N. 


rad 


(C28) 


(C29) 


The solution for the case where there is no radiation or subcooling is given in reference 1 
as 


N fb - 0.41(RaA) 1//4 


When equation (C16) is used, 



(C30) 


A 1 / 4 = 1 _ 0.41 

(RaA) 1 / 4 N fb 


(C31) 


Thus, 



1 - i_ 
0.41 


_3 

F 


1 - f" 


2f 


1_ dlA 

F 2 */ 


1/4 


(C32) 


Changing the form of equation (C26) by expressing B/A 2 ^ 4 in terms of Nusselt numbers 
through the use of equations (C16) and (C17) yields 
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B 

A 3 / 4 


N rad - N tcl 9 
RaA 


(RaA) 


r 3/4 


N rad - N tcl 0 


(RaA) 


1/4 


which from equation (C31) reduces to 


(C33) 



0.41 



Thus, substituting equation (C34) into equation (C26) gives 


N rad - N tci e 


N 


fb 


1 

0.41 


4 t 1 - f 2 1 dF 
F 2f p 2 df 


2 \ 3 / 4 

A _ 1 - f 1 dF ) 

F 2f p 2 df J 


(C34) 


(C35) 


A small correction factor must be applied to the conduction term in the total Nusselt num- 
ber. This correction accounts for the triangular patches shown in figure 1 (p. 2). The 
correction factor, which is derived in reference 1, is given by 


(f) ff ( 1 - f 2 ) + 2 yg - n 1 
° 1 ' 2 yjj 1-f 2 


(C36) 


Numerical Procedure 

The important relations of the preceding section are 

N tot = ^fb + N rad < c3r| ) 


where 
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<p = 


N* 


N tot ^rad 


N. 


fb 


N. 


fb 


and is given parametrically in terms of f by 

„(*) = <p = swt 1 - f2 > (- 1 . lj* ± «Y' /4 

0.41 V F 2f p 2 df I 

and where 9 is related parametrically to radiation and subcooling by 


(C38) 


N. 


rad 


N. 


N tcl 0 1 

—EL = ^(f) = 1 


A , 1 - f 1 dF 


2f 


-,2 df 


fb 


0.41 


2 \ ^/4 

1 - r 1 dF 


_3 

F 


2f 


,2 df 


(C35) 


The function F is given by equation (C8), and £ by equation (C36). These equations 
constitute the general solution of the problem, since a method for computing <p when 
N ra d> Njjy N and 9 are given is apparent. Simply stated, cp is a unique function of 
the grouping (N ra£ j - The numerical procedure is to substitute values of f in 

the range 0 to 1 into the right sides of equations (C38) and (C35). In this way, the com- 
plete functional dependence of (p on radiation and subcooling can be determined. The 
rather involved algebraic manipulations of the preceding section have thereby eliminated 
any need for numerical iteration. The calculations are presented graphically in figure 3 

TABLE in. - TABULATED NUMERICAL RESULTS SHOWING EFFECT 
OF SUBCOOLING AND RADIATION ON DOME SPACING AND 
OVERALL HEAT-TRANSFER COEFFICIENT 


Ratio of vapor -dome 
radius to cell radius, 

f = VR 0 

0.40 
.45 
. 50 
. 55 
.60 
.65 
. 673 + 

.700 

.750 

.790 


Individual heat-transfer 
coefficient group, 

^ = ( h rad “ h tcl^ h fb 

- 1.266 
- 1 . 235 
- 1 . 165 
- 1.032 
-. 7926 
-.3431 
.0000 
.6098 
3.471 
41.56 


Overall heat-transfer 
coefficient group, 

^ = ( h tot ' h rad^ h fb 

1 . 262 
1.252 
1.228 
1 . 188 
1 . 130 
1.049 
1.000 
. 9355 
. 7598 
.3877 
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(p. 13) and selected values are tabulated in table in. 

When values of f = Rj/Rg less than 0. 4 are substituted into the parametric equa- 
tions, a new branch solution is obtained (shown by the dotted line in fig. 3). This branch 
is discussed in the RESULTS section; it is the degenerate case of pure free convection. 


Also, for a value of f = 0. 793 the gap thickness becomes zero and the value of 
(h ra d - htci$)/hfb becomes infinite. At the same time, the ordinate (h^ - h ra( j)/hfk 


becomes zero. 


Thus, for infinite radiation, h^ = h ra( j, 


which is physically correct. 


However, why infinite radiation should correspond to f = 0. 793 is unresolved. 
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APPENDIX D 


SAMPLE CALCULATION 


To illustrate the heat-transfer relations obtained, a numerical example of the film 


boiling of liquid nitrogen is considered. 


•Temperature of 
bulk liquid, 

T|> 

114° R 


-Saturation temper- 
ature of liquid, 

T sat» 

140° R 

Hot wall Wall temperature, 

V 

2740° R 

Figure 5. - Film boiling of liquid nitrogen. 



The various temperatures are shown in figure 5. 
These temperatures were chosen because the 
saturation temperature of nitrogen at atmos- 
pheric pressure is very nearly 140° R and the 
triple point is 114° R. A wall temperature of 
2740° R was selected arbitrarily to assess the 
effect of radiation. To evaluate the heat- 
transfer coefficients the vapor properties are 
evaluated at the film temperature: 


T f = 


T + T 
w s 


2740 + 140 _ 144 qO r _ g8Q o F 


Physical Properties 


The following vapor properties were taken from Eckert and Drake (ref. 13, p. 506) 
and converted to consistent units: 


p = 0. 0267 lb mass/ft 3 
Cp = 0. 2681 Btu/lb mass 
p = 23. 41xl0 -6 lb mass /(sec) (ft) 
k = 0. 09X10" 4 Btu/(sec)(ft)(°R) 

The liquid properties used were taken from reference 14 and evaluated at T . The 

s 

number of the section in which it appears is given after each property. 

Pj = 50. 4 lb mass/ft 3 (Sec. 1. 004) 
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= 1.055x10“^ lb mass/(sec)(ft) 

kj = 2. 24x10“ 5 Btu/(sec)(ft)(°R) 

C } = 0. 49 Btu/(lb mass)(°R) 

P, l 

X = 85. 81 Btu/lb mass 
ct = 6. 03x10 ^ lb force/ft 
/3^ = 3. 23xl0~ 3 (°R) (calculated) 


(Sec. 10. 004) 
(Sec. 3.004) 
(Sec. 4. 004) 
(Sec. 5.004) 
(Sec. 9.004) 
(Sec. 1.004) 


Evaluation of Saturated Film- Boiling Coefficient 


The saturated film-boiling coefficient is evaluated as 


h ffc) = 0.41 


k 3 A*pg(p z - p) 


P-(T, 


w 


T s )i 


1/4 


where 


l = 



6. 03x10 


-4 


1/2 


50.4 - 0.0267 


= 3. 46x10 ft 


and 


A'* = A f 1 + — -P 
20 


C (T -TP 

" V W S' 


= 85. 81 fl + — 0. 2681 X 2600^ = 74g< Q1 Btu/lb 


20 


85. 81 


Therefore, 


h fb = 0.41 


(p. 09x10 4 ) 3 x 748. 01 x 0. 0267 x 32. 2 x (50. 4 - 0. 0267) 


“ 1 1/4 


23. 41X10" 6 x 2600 x 3. 46xl0~ 3 


x 3600 
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or 


h fb = 27. 00 Btu/ (hr) (ft 2 ) (°R) 

Evaluation of Radiation Coefficient 

The plate and liquid emissivities are assumed to be 1 in the following evaluation of 
h rad : 


rp4 _ rp4 

h rad = °* 1713x l 0 " 8 — = 0. 1713 

T w “ T s 


27. 4 4 - 1. 4 4 
2600 


= 37. 50 Btu / (hr) (ft 2 ) (°R) 


Subcooling Parameter 

For the film boiling of liquid nitrogen, the subcooling parameter is evaluated as 

e = = o. oi 

T w - T g 2600 


Liquid Free-Convection Coefficient 


One of the major questions left unanswered by the present analysis is what to use for 

the turbulent liquid free-convection coefficient, since the vapor bubbles induce strong 

turbulence in the liquid. Experimental data in this area would be most helpful. For the 

purposes of this calculation, a coefficient based on ordinary turbulent free convection 

from a horizontal plate given by McAdams (ref. 15, p. 180) is used. The liquid properties 

were evaluated at T , though McAdams recommends using the liquid film temperature: 

s 


h tcl = °- 14 


k i p ; gl V T s - T b> 

2 

H 


1/3 


Pr, 


The Prandtl number is given by 
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C p, l H = 0.49 x 1. 05 5xlQ~ 4 = 2 31 
k Z 2. 24x10" 5 


Therefore, 


h tcl " °* 14 


(2. 24x10" 5 f x 50. 4 2 x 32. 2 x 3^ 23x10" 3 >< 26 x 2. 31 


11/3 


x 3600 


( 1 . 055X10 -4 ) 


= 127. 15 B tu/(hr ) (ft 2 ) (°R) 


Evaluation of Overall Heat -Transfer Coefficient 


To evaluate the overall heat-transfer coefficient 


h rad ~ h tcl g = 37.50 - 127. 15 x 0.01 = ^ 342 
11^ 27. 00 

is first computed. From figure 3 (p. 13), the value of cp is 0.871. Thus, 

h tot = ^ h fb + h rad = °* 871 X 27 ' 00 + 37 * 50 = 61 * 02 Btu/(hr)(ft 2 )(°R) 
and from equation (31) 

q tot = 61 * 02 x 2600 = 1 - 59xl ° 5 Btu/(hr)(ft 2 ) 

The simple formula given by equation (38) predicts 

btot = h fb + °- 88h rad + °* 12h tcl 0 = 27 ‘ 00 + °* 88 x 37 - 5 + 0. 12 x 127. 15 x 0. 01 
= 60. 16 Btu/(hr ) (ft 2 ) (°R) 


This value is very close to the graphical solution, indicating the utility of the simple 
formula. 
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APPENDIX E 

TRANSFORMATION OF VARIABLES 

The optimization constraints, equations (22) and (23), are rewritten in terms of the 
parameter f for convenience in the mathematical operations, as given by equation (Cl). 
The total heat -transfer coefficient can be rewritten in terms of the parameter f as 

h tot^ R 0’ R l) “ h tot^ R 0’ ( El ) 

The functional forms of h^ on the left and right side of this equation are, of course, 
different. 

The optimal constraint, equation (23), can now be rewritten as 

^k.®rHo)] Ro ^|-M.Ho)] Ro = <- < E2 > 


or 




(E3) 


which is used in arriving at equation (C14). 

Now, by the chain rule of differentiation, the optimal constraint, equation (22), can 
be rewritten as 


a 

9R o 



dR o 




9f 

R o aR o 


= 0 


(E4) 


However, because of equation (E3), equation (E4) simplifies to 


9 

aR o 




(E5) 


which is equation (Cl). 
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